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Abstract 

C/3 . 

, We represent planning as a set of loosely coupled network flow problems, where each 

network corresponds to one of the state variables in the planning domain. The network 
nodes correspond to the state variable values and the network arcs correspond to the value 
^ \ transitions. The planning problem is to find a path (a sequence of actions) in each network 

such that, when merged, they constitute a feasible plan. In this paper we present a num- 
ber of integer programming formulations that model these loosely coupled networks with 
' varying degrees of flexibility. Since merging may introduce exponentially many ordering 

constraints we implement a so-called branch-and-cut algorithm, in which these constraints 
are dynamically generated and added to the formulation when needed. Our results are very 
promising, they improve upon previous planning as integer programming approaches and 
lay the foundation for integer programming approaches for cost optimal planning. 



in 

o 



^ ■ 1. Introduction 



While integer programming^ approaches for automated planning have not been able to 
scale well against other compilation approaches (i.e. satisfiability and constraint satisfac- 
tion), they have been extremely successful in the solution of many real- world large scale 
optimization problems. Given that the integer programming framework has the potential 
to incorporate several important aspects of real-world automated planning problems (for 
example, numeric quantities and objective functions involving costs and utilities), there is 
significant motivation to investigate more effective integer programming formulations for 
classical planning as they could lay the groundwork for large scale optimization (in terms 
of cost and resources) in automated planning. In this paper, we study a novel decomposi- 
tion based approach for automated planning that yields very effective integer programming 
formulations. 



1. We use the term integer programming to refer to integer linear programming unless stated otherwise. 
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Decomposition is a general approach to solving problems more efRcicntly. It involves 
breaking a problem up into several smaller subproblems and solving each of the subprob- 
lems separately. In this paper we use decomposition to break up a planning problem into 
several interacting (i.e. loosely coupled) components. In such a decomposition, the planning 
problem involves both finding solutions to the individual components and trying to merge 
them into a feasible plan. This general approach, however, prompts the following questions: 
(1) what are the components, (2) what are the component solutions, and (3) how hard is it 
to merge the individual component solutions into a feasible plan? 

1.1 The Components 

We let the components represent the state variables of the planning problem. Figure 1 
illustrates this idea using a small logistics example, with one truck and a package that 
needs to be moved from location 1 to location 2. There are a total of five components in 
this example, one for each state variable. We represent the components by an appropriately 
defined network, where the network nodes correspond to the values of the state variable 
(for atoms this is T = true and F = false), and the network arcs correspond to the value 
transitions. The source node in each network, represented by a small in-arc, corresponds 
to the initial value of the state variable. The sink node(s), represented by double circles, 
correspond to the goal value(s) of the state variable. Note that the effects of an action 
can trigger value transitions in the state variables. For example, loading the package at 
location 1 makes the atom pack-in-truck true and pack-at-locl false. In addition, loading 
the package at location 1 requires that the atom truck-at-locl is true. 

While the idea of components representing the state variables of the planning problem 
can be used with any state variable representation, it is particularly synergistic with multi- 
vahicd state variables. MTilti-valucd state variables provide a more compact representation 
of the planning problem than their binary-valued counterparts. Therefore, by making the 
conversion to multi-valued state variables we can reduce the number of components and 
create a better partitioning of the constraints. Figure 2 illustrates the use of multi-valued 
state variables on our small logistics example. There are two multi-valued state variables 
in this problem, one to characterize the location of the truck and one to characterize the 
location of the package. In our network representation, the nodes correspond to the state 
variable values (1 = at-locl, 2 = at-loc2, and t = in-truck), and the arcs correspond to the 
value transitions. 

1.2 The Component Solutions 

We let the component solutions represent a path of value transitions in the state variables. 
In the networks, nodes and arcs appear in layers. Each layer represents a plan period 
in which, depending on the structure of the network, one or more value transitions can 
occur. The networks in Figures 1 and 2 each have three layers (i.e. plan periods) and their 
structure allows values to persist or change exactly once per period. The layers are used to 
solve the planning problem incrementally. That is, we start with one layer in each network 
and try to solve the planning problem. If no plan is found, all networks are extended by 
one extra layer and a new attempt is made to solve the planning problem. This process is 
repeated until a plan is found or a time limit is reached. In Figures 1 and 2, a path (i.e. 
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loci loc2 




Figure 1: Logistics example broken up into five components (binary- valued state variables) 
that are represented by network flow problems. 




Figure 2: Logistics example broken up into two components (multi-valued state variables) 
that are represented by network flow problems. 



a solution) from the source node to one of the sink nodes is highlighted in each network. 
Since the execution of an action triggers value transitions in the state variables, each path 
in a network corresponds to a sequence of actions. Consequently, the planning problem 
can he thought of as a collection of network flow problems where the problem is to find a 
path (i.e. a sequence of actions) in each of the networks. However, interactions between 
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the networks impose side constraints on the network flow problems, which complicate the 
solution process. 

1.3 The Merging Process 

We solve these loosely coupled networks using integer programming formulations. One 

design choice wc make is that we expand all networks (i.e. components) together, so the 
cost of finding solutions for the individual networks as well as merging them depends on the 
difficulty of solving the integer programming formulation. This, in turn, typically depends 
on the size of the integer programming formulation, which is partly determined by the 
number of layers in each of the networks. The simplest idea is to have the number of 
layers of the networks equal the length of the plan, just as in sequential planning where the 
plan length equals the number of actions in the plan. In this case, there will be as many 
transitions in the networks as there are actions in the plan, with the only difference that a 
sequence of actions corresponding to a path in a network could contain no-op actions. 

An idea to reduce the required number of layers is by allowing multiple actions to be 
executed in the same plan period. This is exactly what is done in Graphplan (Blum & Furst, 
1995) and in other planners that have adopted the Graphplan-style definition of parallelism. 
That is, two actions can be executed in parallel (i.e. in the same plan period) as long as 
they are non-interfering. In our formulations we adopt more general notions of parallelism. 
In particular, wc relax the strict relation between the mimbcr of layers in the networks 
and the length of the plan by changing the network representation of the state variables. 
For example, by allowing multiple transitions in each network per plan period we permit 
interfering actions to be executed in the same plan period. This, however, raises issues 
about how solutions to the individual networks arc searched and how they are combined. 
When the network representations for the state variables allow multiple transitions in each 
network per plan period, and thus become more flexible, it becomes harder to merge the 
solutions into a feasible plan. Therefore, to evaluate the tradeoffs in allowing such flexible 
representations, we look at a variety of integer programming formulations. 

We refer to the integer programming formulation that uses the network representation 
shown in Figures 1 and 2 as the one state ehange model, because it allows at most one 
transition (i.e. state change) per plan period in each state variable. Note that in this network 
representation a plan period mimics the Graphplan-style parallelism. That is, two actions 
can be executed in the same plan period if one action does not delete the precondition or 
add-effect of the other action. A more flexible representation in which values can change at 
most once and persist before and after each change we refer to as the generalized one state 
change model. Clearly, we can increase the number of changes that we allow in each plan 
period. The representations in which values can change at most twice or k times, we refer 
to as the generalized two state change and the generalized k state change model respectively. 
One disadvantage with the generalized k state change model is that it creates one variable 
for each way to do k value changes, and thus introduces exponentially many variables per 
plan period. Therefore, another network representation that we consider allows a path of 
value transitions in which each value can be visited at most once per plan period. This 
way, we can limit the number of variables, but may introduce cycles in our networks. The 
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integer programming formulation that uses this representation is referred to as the state 
change path model. 

In general, by allowing multiple transitions in each network per plan period (i.e. layer), 
the more complex the merging process becomes. In particular, the merging process checks 
whether the actions in the solutions of the individual networks can be linearized into a 
feasible plan. In our integer programming formulations, ordering constraints ensure feasible 
linearizations. There may, however, be exponentially many ordering constraints when we 
generalize the Graphplan-style parallelism. Rather than inserting all these constraints in the 
integer programming formulation up front, we add them as needed using a branch-and-cut 
algorithm. A branch-and-cut algorithm is a branch-and-bound algorithm in which certain 
constraints are generated dynamically throughout the branch-and-bound tree. 

We show that the performance of our integer programming (IP) formulations show new 
potential and are competitive with SATPLAN04 (Kautz, 2004). This is a significant result 
because it forms a basis for other more sophisticated IP-based planning systems capable of 
handling numeric constraints and non-uniform action costs. In particular, the new potential 
of our IP formulations has led to their successful use in solving partial satisfaction planning 
problems (Do, Benton, van den Briel, &; Kambhampati, 2007). Moreover, it has initiated a 
new line of work in which integer and linear programming are used in heuristic state-space 
search for automated planning (Benton, van den Briel, &: Kambhampati, 2007; van den 
Briel, Benton, Kambhampati, & Vossen, 2007). 

The remainder of this paper is organized as follows. In Section 2 we provide a brief 
background on integer programming and discuss some approaches that have used integer 
programming to solve planning problems. In Section 3 we present a series of integer pro- 
gramming formulations that each adopt a different network representation. We describe how 
we set up these loosely coupled networks, provide the corresponding integer programming 
formulation, and discuss the different variables and constraints. In Section 4 we describe 
the branch-and-cut algorithm that is used for solving these formulations. We provide a 
general background on the branch-and-cut concept and show how we apply it to our for- 
mulations by means of an example. Section 5 provides experimental results to determine 
which characteristics in our approach have the greatest impact on performance. Related 
work is discussed in Section 6 and some conclusions are given in Section 7. 

2. Background 

Since our formulations are based on integer programming, we briefly review this technique 
and discuss its use in planning. A mixed integer program is represented by a linear objective 
function and a set of linear inequalities: 

minjcx : Ax > b, xi, Xp >Q and integer, Xp+i, x„ > 0}, 

where A is an (m x n) matrix, c is an n-dimensional row vector, b is an m-dimensional 
column vector, and x an n-dimensional column vector of variables. If all variables are 
continuous {p = 0) we have a linear program, if all variables are integer (p = n) we have 
an integer program, and if xi, ...,Xp € {0, 1} we have a mixed 0-1 program. The set S = 
{xi,...,Xp > and integer, Xp+i, x„ > : Ax > 6} is called the feasible region, and an 
n-dimensional column vector x is called a feasible solution if x € 5. Moreover, the function 
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cx is called the objective function, and the feasible solution x* is called an optimal solution 
if the objective function is as small as possible, that is, cx* = minjca; : x € S} 

Mixed integer programming provides a rich modeling formalism that is more general than 
propositional logic. Any propositional clause can be represented by one linear inequality in 
0-1 variables, but a single linear inequality in 0-1 variables may require exponentially many 
clauses (Hooker, 1988). 

The most widely used method for solving (mixed) integer programs is by applying a 
branch-and-bound algorithm to the linear programming relaxation, which is much easier 
to solve^. The linear programming (LP) relaxation is a linear program obtained from the 
original (mixed) integer program by relaxing the integrality constraints: 

min{cx : Ax > b, xi, x„ > 0} 

Generally, the LP relaxation is solved at every node in the branch-and-bound tree, until 
(1) the LP relaxation gives an integer solution, (2) the LP relaxation value is inferior to the 
current best feasible solution, or (3) the LP relaxation is infeasible, which implies that the 
corresponding (mixed) integer program is infeasible. 

An ideal formulation of an integer program is one for which the solution of the lin- 
ear programming relaxation is integral. Even though every integer program has an ideal 
formulation (Wolsey, 1998), in practice it is very hard to characterize the ideal formula- 
tion as it may require an exponential number of inequalities. In problems where the ideal 
formulation cannot be determined, it is often desirable to find a strong formulation of 
the integer program. Suppose that the feasible regions Pi = {x £ i?" : A^x > bi} and 
P2 = {x e RT' : A2X > 62} describe the linear programming relaxations of two IP formula- 
tions of a problem. Then we say that formulation for Pi is stronger than formulation for P2 
if Pi d P2- That is, the feasible region Pi is subsumed by the feasible region P2. In other 
words Pi improves the quality of the linear relaxation of P2 by removing fractional extreme 
points. 

There exist numerous powerful software packages that solve mixed integer programs. In 
our experiments we make use of the commercial solver CPLEX 10.0 (Inc., 2002), which is 
currently one of the best LP/IP solvers. 

The use of integer programming techniques to solve artificial intelligence planning prob- 
lems has an intuitive appeal, especially given the success IP has had in solving similar types 
of problems. For example, IP has been used extensively for solving problems in transporta- 
tion, logistics, and manufacturing. Examples include crew scheduling, vehicle routing, and 
production planning problems (Johnson, Nemhauser, &; Savelsbergh, 2000). One potential 
advantage is that IP techniques can provide a natural way to incorporate several important 
aspects of real- world planning problems, such as numeric constraints and objective functions 
involving costs and utilities. 

Planning as integer programming has, nevertheless, received only limited attention. One 
of the first approaches is described by Bylander (1997), who proposes an LP heuristic for 
partial order planning algorithms. While the LP heuristic helps to reduce the number of 
expanded nodes, the evaluation is rather time-consuming. In general, the performance of 

2. While the integer programming problem is A'^P-complete (Garey & Johnson, 1979) the linear program- 
ming problem is polynomially solvable (Karmarkar, 1984). 
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IP often depends on the structure of the problem and on how the problem is formulated. 
The importance of developing strong IP formulations is discussed by Vossen ct al. (1999), 
who compare two formulations for classical planning: (1) a straightforward formulation 
based on the conversion of the propositional representation by SATPLAN which yields 
only mediocre results, and (2) a less intuitive formulation based on the representation of 
state transitions which leads to considerable performance improvements. Several ideas that 
further improve formulation based on the representation of state transitions are described by 
Dimopoulos (2001). Some of these ideas are implemented in the IP-based planner Optiplan 
(van den Briel k, Kambhampati, 2005). Approaches that rely on domain-specific knowledge 
are proposed by Bockmayr and Dimopoulos (1998, 1999). By exploiting the structure of 
the planning problem these IP formulations often provide encouraging results. The use of 
LP and IP has also been explored for non-classical planning. Dimopoulos and Gerevini 
(2002) describe an IP formulation for temporal planning and Wolfman and Weld (1999) 
use LP formulations in combination with a satisfiability-based planner to solve resource 
planning problems. Kautz and Walser (1999) also solve resource planning problems, but 
use domain-specific IP formulations. 

3. Formulations 

This section describes four IP formulations that model the planning problem as a collection 
of loosely coupled network flow problems. Each network represents a state variable, in which 
the nodes correspond to the state variable values, and the arcs correspond to the value 
transitions. The state variables are based on the SAS-I- planning formalism (Backstrom & 
Nebel, 1995), which is a planning formalism that uses multi- valued state variables instead 
of binary-valued atoms. An action in SAS-I- is modeled by its pre-, post- and prcvail- 
conditions. The pre- and post-conditions express which state variables arc changed and 
what values they must have before and after the execution of the action, and the prevail- 
conditions specify which of the unchanged variables must have some specific value before 
and during the execution of an action. A SAS-I- planning problem is described by a tuple 
n = (C, A, so,s^) where: 

• C = {ci, c„} is a finite set of state variables, where each state variable c G C has an 

associated domain Vc and an implicitly defined extended domain = V^U{n}, where 
u denotes the undefined value. For each state variable c € C, s[c] denotes the value of 
c in state s. The value of c is said to be defined in state s if and only if s[c] / u. The 
total state space S = Vc^ x ... x Vc„ and the partial state space 5+ = x ... x 
are implicitly defined. 

• ^ is a finite set of actions of the form {pre, post, prev), where pre denotes the pre- 
conditions, post denotes the post-conditions, and prev denotes the prevail-conditions. 
For each action a G A, pre[c],post[c] and prev[c] denotes the respective conditions on 
state variable c. The following two restrictions are imposed on all actions: (1) Once 
the value of a state variable is defined, it can never become undefined. Hence, for all 
c G C, if pre[c] ^ u then pre[c] ^ post[c] ^ u; (2) A prevail- and post-condition of 
an action can never define a value on the same state variable. Hence, for all c G C, 
either post[c] = u or prev[c] = u oi both. 
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• sq e S denotes the initial state and G 5"*" denotes the goal state. While SAS+ 
planning allows the initial state and goal state to be both partial states, we assume 
that So is a total state and is a partial state. We say that state s is satisfied by 
state t if and only if for all c G C we have s[c] = u or s[c] = t[c]. This implies that if 
s* [c] = u for state variable c, then any defined value / G K satisfies the goal for c. 

To obtain a SAS+ description of the planning problem we use the translator component 
of the Fast Downward planner (Helmert, 2006). The translator is a stand-alone component 
that contains a general purpose algorithm which transforms a propositional description 
of the planning problem into a SAS+ description. The algorithm provides an efficient 
grounding that minimizes the state description length and is based on the preprocessing 
algorithm of the MIPS planner (Edelkamp & Helmert, 1999). 

In the remainder of this section we introduce some notation and describe our IP formu- 
lations. The formulations are presented in such a way that they progressively generalize the 
Graphplan-style parallelism through the incorporation of more flexible network representa- 
tions. For each formulation we will describe the underlying network, and define the variables 
and constraints. We will not concentrate on the objective function as much because the 
constraints will tolerate only feasible plans. 

3.1 Notation 

For the formulations that are described in this paper we assume that the following infor- 
mation is given: 

• C: a set of state variables; 

• Vc'. a set of possible values (i.e. domain) for each state variable c G C; 

• Ec- & set of possible value transitions for each state variable c G C; 

• Gc = {Vc, Ec) : a directed domain transition graph for every c G C; 

State variables can be represented by a domain transition graph, where the nodes correspond 
to the possible values, and the arcs correspond to the possible value transitions. An example 
of the domain transition graph of a variable is given in Figure 3. While the example depicts 
a complete graph, a domain transition graph does not need to be a complete graph. 
Furthermore, we assume as given: 

• E^ C Ec represents the effect of action a in c; 

• C Vc represents the prevail condition of action a in c; 

• Ac := {a £ A : \Ec \ > 0} represents the actions that have an effect in c, and Af{e) 
represents the actions that have the effect e in c; 

• AY := {a G a : \Vc\ > 0} represents the actions that have a prevail condition in c, 

and A]^{f) represents the actions that have the prevail condition / in c; 

• C"" := {c e C : a e Ac U A^} represents the state variables on which action a has an 
effect or a prevail condition. 
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Figure 3: An example of a domain transition graph, where Vc = {f,g, h} are the possible 
values (states) of c and Ec = {{f,g),{f,h),{g,f),{g,h),{h,f),{h,g)} are the 
possible value transitions in c. 



Hence, each action is defined by its effects (i.e. pre- and post-conditions) and its prevail 
conditions. In SAS-I- planning, actions can have at most one effect or prevail condition in 
each state variable. In other words, for each a e A and c G C, we have that and are 
empty or -|- |V^"| < 1. An example of how the effects and prevail conditions affect one 
or more domain transition graphs is given in Figure 4. 






Figure 4: An example of how action effects and prevail conditions are represented in a 

domain transition graph. Action a has implications on three state variables = 
{ci, C2, C3}. The effects of a are represented by E^^ = {(/, g)} and E^^ = {{h,f )}, 
and the prevail condition of a is represented by V^" = {/i}. 



In addition, we use the following notation: 

• denote the in-arcs of node / in the domain transition graph Gc, 

• V~{f): to denote the out-arcs of node / in the domain transition graph Gc, 

• ^ckif)'- denote paths of length k in the domain transition graph Gc that end at 



node /. Note that P+i(/) = V+{f). 



• P^kif)'- t° denote paths of length k in the domain transition graph Gc that start at 
node /. Note that Pci{f) = V-{f). 
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• -^cfe(/)- t° denote paths of length k in the domain transition graph Gc that visit node 
/, but that do not start or end at /. 

3.2 One State Change (ISC) Formulation 

Our first IP formulation incorporates the network representation that we have seen in 
Figures 1 and 2. The name one state change relates to the number of transitions that we 
allow in each state variable per plan period. The restriction of allowing only one value 
transition in each network also restricts which actions we can execute in the same plan 
period. It happens to be the case that the network representation of the ISC formulation 
incorporates the standard notion of action parallelism which is used in Graphplan (Blum 
&: Furst, 1995). The idea is that actions can be executed in the same plan period as long 
as they do not delete the precondition or add-effcct of another action. In terms of value 
transitions in state variables, this is saying that actions can be executed in the same plan 
period as long as they do not change the same state variable (i.e. there is only one value 
change or value persistence in each state variable). 

3.2.1 State Change Network 

Figure 5 shows a single layer (i.e. period) of the network which underlies the ISC formula- 
tion. If we set up the IP formulation with T plan periods, then there will be T + 1 layers 
of nodes and T layers of arcs in the network (the zeroth layer of nodes is for the initial 
state and the remaining T layers of nodes and arcs are for the successive plan periods). For 
each possible state transition there is an arc in the state change network. The horizontal 
arcs correspond to the persistence of a value, and the diagonal arcs correspond to the value 
changes. A solution path to an individual network follows the arcs whose transitions are 
supported by the action effect and prevail conditions that appear in the solution plan. 

ISC network 




Period t 



Figure 5: One state change (ISC) network. 

3.2.2 Variables 

We have two types of variables in this formulation: action variables to represent the execu- 
tion of an action, and arc flow variables to represent the state transitions in each network. 
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We use separate variables for changes in a state variable (the diagonal arcs in the ISC 
network) and for the persistence of a value in a state variable (the horizontal arcs in the 
ISC network). The variables are defined as follows: 

• e {0, 1}, ioT a e A, 1 < t < T; xf is equal to 1 if action a is executed at plan period 
t, and otherwise. 

• yc,f,t € {0, 1}, for c e C , f £ Vc, I < t < T; yc,f,t is equal to 1 if the value / of state 
variable c persists at period t, and otherwise. 

• yc,e,t G {0, 1}, for c G C, e G Ec, I < t < T; yc,e,t is equal to 1 if the transition e e Ec 
in state variable c is executed at period t, and otherwise. 



3.2.3 Constraints 

There are two classes of constraints. We have constraints for the network flows in each 
state variable network and constraints for the action effects that determine the interactions 
between these networks. The ISC integer programming formulation is: 

• State change flows for all c G C, / G K 

Ef 1 if / = sofcl . 
y^'^'^ + y^J'^ = [ otherwise. 

yc,e,t+l + ycj,t+l = J2 2/c,e,t + yc,f,t for 1 < t < T - 1 (2) 

eevrif) eev+{f) 

Y yc,e,T + yc,f,T = 1 if/ = S*[c] (3) 

eev+{f) 

• Action implications for all c G C, 1 < i < T 

J2 < = yc,e,t ioreeEc (4) 
x^ < y,j,t ioraeAJeV^ (5) 



Constraints (1), (2), and (3) are the network flow constraints for state variable c G C. 
Constraint (1) ensures that the path of state transitions begins in the initial state of the 
state variable and constraint (3) ensures that, if a goal exists, the path ends in the goal 
state of the state variable. Note that, if the goal value for state variable c is undefined 
(i.e. s<,[c] = u) then the path of state transitions may end in any of the values / G V^. 
Hence, we do not need a goal constraint for the state variables whose goal states s*[c] are 
undefined. Constraint (2) is the flow conservation equation and enforces the continuity of 
the constructed path. 

Actions may introduce interactions between the state variables. For instance, the effects 
of the load action in our logistics example affect two different state variables. Actions link 
state variables to each other and these interactions are represented by the action implication 
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constraints. For each transition e G Ec, constraints (4) link the action execution variables 
that have e as an effect (i.e. e € E^) to the arc flow variables. For example, if an action 

with effect e G is executed, then the path in state variable c must follow the arc 
represented by yc,e,t- Likewise, if we choose to follow the arc represented by yc,e,t, then 
exactly one action with e e E^ must be executed. The summation on the left hand side 
prevents two or more actions from interfering with each other, hence only one action may 
cause the state change e in state variable c at period t. 

Prevail conditions of an action link state variables in a similar way as the action effects 
do. Specifically, constraint (5) states that if action a is executed at period t {xf = 1), then 
the prevail condition / G V^" is required in state variable c at period t {ycj,t = 1)- 

3.3 Generalized One State Change (GISC) Formulation 

In our second formulation we incorporate the same network representation as in the ISC 
formulation, but adopt a more general interpretation of the value transitions, which leads 
to an unconventional notion of action parallelism. For the GISC formulation we relax the 
condition that parallel actions can be arranged in any order by requiring a weaker condition. 
We allow actions to be executed in the same plan period as long as there exists some ordering 
that is feasible. More specifically, within a plan period a set of actions is feasible if (1) there 
exists an ordering of the actions such that all preconditions are satisfied, and (2) there is 
at most one state change in each of the state variables. This generalization of conditions 
is similar to what Rintanen, Heljanko and Niemela (2006) refer to as the 3-step semantics 
semantics. 

To illustrate the basic concept, let us again examine our small logistics example intro- 
duced in Figure 1. The solution to this problem is to load the package at location 1, drive the 
truck from location 1 to location 2, and unload the package at location 2. Clearly, this plan 
would require three plan periods under Graphplan-style parallelism as these three actions 
interfere with each other. If, however, we allow the load at loci and the drive loci — )■ loc2 
action to be executed in the same plan period, then there exists some ordering between 
these two actions that is feasible, namely load the package at the location 1 before driving 
the truck to location 2. The key idea behind this example should be clear: while it may 
not be possible to find a set of actions that can be linearized in any order, there may nev- 
ertheless exist some ordering of the actions that is viable. The question is, of course, how 
to incorporate this idea into an IP formulation. 



Load at loci 




Figure 6: Logistics example represented by network flow problems with generalized arcs. 
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This example illustrates that we are looking for a set of constraints that allow sets of 
actions for which: (1) all action preconditions are met, (2) there exists an ordering of the 
actions at each plan period that is feasible, and (3) within each state variable, the value is 
changed at most once. The incorporation of these ideas only requires minor modifications 
to the ISC formulation. Specifically, we need to change the action implication constraints 
for the prevail conditions and add a new set of constraints which we call the ordering 
implication constraints. 

3.3.1 State Change Network 

The minor modifications are revealed in the GISC network. While the network itself is 

identical to the ISC network, the interpretation of the transition arcs is somewhat different. 
To incorporate the new set of conditions, we implicitly allow values to persist (the dashed 
horizontal arcs in the CISC network) at the tail and head of each transition arc. The 
interpretation of these implicit arcs is that in each plan period a value may be required as 
a prevail condition, then the value may change, and the new value may also be required as 
a prevail condition as shown in Figure 7. 



GISC network Generalized state change arc 




Period t Period t 



Figure 7: Generalized one state change (GISC) network. 

3.3.2 Variables 

Since the GISC network is similar to the ISC network the same variables are used, thus, 
action variables to represent the execution of an action, and arc flow variables to represent 
the flow through each network. The difference in the interpretation of the state change arcs 
is dealt with in the constraints of the CISC formulation, and therefore does not introduce 
any new variables. For the variable definitions, we refer to Section 3.2.2. 

3.3.3 Constraints 

We now have three classes of constraints, that is, constraints for the network flows in each 
state variable network, constraints for linking the flows with the action effects and prevail 
conditions, and ordering constraints to ensure that the actions in the plan can be linearized 
into a feasible ordering. 
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The network flow constraints for the GISC formulation are identical to those in the ISC 
formulation given by (l)-(3). Moreover, the constraints that link the flows with the action 
effects are equal to the action effect constraints in the ISC formulation given by (4). The 
CISC formulation differs from the ISC formulation in that it relaxes the condition that 
parallel actions can be arranged in any order by requiring a weaker condition. This weaker 
condition affects the constraints that link the flows with the action prevail conditions, and 
introduces a new set of ordering constraints. These constraints of the CISC formulation 
are given as follows: 

• Action implications for all cGC, l<t<T 

xt < yc,f,t+ y'^'^'t^ H 2/c,e,t ioraeAJeV^^ (6) 

eeVc+(/) eeVc-{f) 

• Ordering implications 

J2 < |'t^(A)| -1 for all cycles AG 0^''"'='= (7) 

aeV(A) 

Constraint (6) incorporates this new set of conditions for which actions can be executed 
in the same plan period. In particular, we need to ensure that for each state variable c, the 
value f & Vc holds if it is required by the prevail condition of action a at plan period t. There 
are three possibilities: (1) The value / holds for c throughout the period. (2) The value / 
holds initially for c, but the value is changed to a value other than / by another action. (3) 
The value / does not hold initially for c, but the value is changed to / by another action. 
In either of the three cases the value / holds at some point in period t so that the prevail 
condition for action a can be satisfied. In words, the value / may prevail implicitly as long 
as there is a state change that includes /. As before, the prevail implication constraints 
link the action prevail conditions to the corresponding network arcs. 

The action implication constraints ensure that the preconditions of the actions in the 
plan are satisfied. This, however, does not guarantee that the actions can be linearized 
into a feasible order. Figure 7 indicates that there are implied orderings between actions. 
Actions that require the value / as a prevail condition must be executed before the action 
that changes / into g. Likewise, an action that changes / into g must be executed before 
actions that require the value ^ as a prevail condition. The state change flow and action 
implication constraints outlined above indicate that there is an ordering between the actions, 
but this ordering could be cyclic and therefore infeasible. To make sure that an ordering 
is acyclic we start by creating a directed implied precedence graph G^'^^'^ = {yP'^<^<^^E^^'^). 
In this graph the nodes a G V'^^^'^ correspond to the actions, that is, V'^^^^ = A, and we 
create a directed arc (i.e. an ordering) between two nodes (a, h) G E^'^^'^ if action a has to be 
executed before action h in time period t, or if h has to be executed after a. In particular, 
we have 

S^'-^ = U (a,6)U U (a, 6) 

{a,b)eAy.A,c&CJeV^,eeE^^: (a,6)eAx A,cGC,geyi',eeE°: 
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The implied ordcrings become immediately clear from Figure 8. The figure on the left 
depicts the first set of orderings in the expression of E^'^^'^. It says that the ordering between 
two actions a and h that are executed in the same plan period is implied if action a requires 
a value to prevail that action h deletes. Similarly, the figure on the right depicts second set 
of orderings in the expression of E^^^'^. That is, an ordering is implied if action a adds the 
prevail condition of 6. 




Figure 8: Implied orderings for the GISC formulation. 

The ordering implication constraints ensure that the actions in the final solution can be 
linearized. They basically involve putting an n-ary mutex relation between the actions that 
arc involved in each cycle. Unfortunately, the number of ordering implication constraints 
grows exponentially in the number of actions. As a result, it will be impossible to solve 
the resulting formulation using standard approaches. We address this complication by 
implementing a branch-and-cut approach in which the ordering implication constraints are 
added dynamically to the formulation. This approach is discussed in Section 4. 

3.4 Generalized k State Change (G/cSC) Formulation 

In the GISC formulation actions can be executed in the same plan period if (1) there exists 
an ordering of the actions such that all preconditions are satisfied, and (2) there occurs at 
most one value change in each of the state variables. One obvious generalization of this 
would be to relax the second condition and allow at most value changes in each state 
variable c, where kc < \Vc\ — I- By allowing multiple value changes in a state variable per 
plan period we, in fact, permit a series of value changes. Specifically, the GfcSC model 
allows series of value changes. 

Obviously, there is a tradeoff between loosening the networks versus the amount of 
work it takes to merge the individual plans. While we have not implemented the GfcSC 
formulation, we provide some insight in this tradeoff by describing and evaluating the GfcSC 
formulation with kc = 2 for all c G C We will refer to this special case as the generalized 
two state change (G2SC) formulation. One reason we restrict ourselves to this special case 
is that the general case of k state changes would introduce exponentially many variables 
in the formulation. There are IP techniques, however, that deal with exponentially many 
variables (Desaulniers, Desrosiers, Sz Solomon, 2005), but we will not discuss them here. 

3.4.1 State Change Network 

The network that underlies the G2SC formulation is equivalent to CISC, but spans an extra 
layer of nodes and arcs. This extra layer allows us to have a series of two transitions per plan 
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period. All transitions are generalized and implicitly allow values to persist just as in the 
GISC network. Figure 9 displays the network corresponding to the G2SC formulation. In 
the G2SC network there are generalized one and two state change arcs. For example, there 
is a generalized one state change arc for the transition {f,g), and there is a generalized two 
state changes arc for the transitions {(/, g), {g, h)}. Since all arcs are generalized, each value 
that is visited can also be persisted. We also allow cyclic transitions, such as, {(/, g), (g, /)} 
if / is not the prevail condition of some action. If we were to allow cyclic transitions in 
which / is a prevail condition of an action, then the action ordering in a plan period can not 
be implied anymore (i.e. the prevail condition on / would either have to occur before the 
value transitions to g, or after it transitions back to /). Thus if there is no prevail condition 
on / then we can safely allow the cyclic transition {(f,g), (g, /)}. 




Figure 9: Generalized two state change (G2SC) network. On the left the subnetwork that 
consists of generalized one state change arcs and no-op arcs, on the right the sub- 
network that consists of the generalized two state change arcs. The subnetwork 
for the two state change arcs may include cyclic transitions, such as, {(/, g), {g, /)} 
as long as / is not the prevail condition of some action. 



3.4.2 Variables 

As before we have variables representing the execution of an action, and variables repre- 
senting the flows over one state change (diagonal arcs) or persistence (horizontal arcs). In 
addition, we have variables representing paths over two consecutive state changes. Hence, 
we have variables for each pair of state changes (/, g, h) such that (/, g) G and {g, h) G Ec- 
We will restrict these paths to visit unique values only, that is, f ^ g, g h, and h ^ f, 
or if / is not a prevail condition of any action then we also allow paths where f = h. The 
variables from the GISC formulation are also used in G2SC formulation. There is, however, 
an additional variable to represent the arcs that allow for two state changes: 

• yc,ei,e2,t G {0,1}, for c € C, (61,62) e Pc,2, I < t < T; yc,ei,e2,t is equal to 1 if 
there exists a value f & Vc and transitions 61,62 € Ec, such that ei G V^^(/) and 
62 G V~{f), in state variable c are executed at period t, and otherwise. 
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3.4.3 Constraints 

We again have our three classes of constraints, which are given as follows: 

• State change flows for all c G C, / G Vc 

yc,e„e„i+ }^ yc,e,i + ycj,i = jo otherwise. 

(ei,e2)eP-2(/) eeVc'if) 

yc,ei,e2,t+l + X] 2/c,e,t+l + yc,f,t+l = 
(ei,e2)eP-2(/) eeVrif) 

Yl yc,ei,e2,t + Yl 2/c,e,t + ycj,t for 1 < t < T - 1 (9) 

(ei,e2)eP+2(/) eeVc+{f) 

J2 yc,ei,e2,T+ yc,e,T + ycJ,T = 1 if{/GS*[c]} (10) 

(ei,e2)eP+2(/) eeV+{f) 

• Action implications for all c G C, 1 <t <T 

Y = yc,e,t + Y yc,ei,e2,t for e G (H) 

aeA-.eeEg (ei,e2)ePc,2:ei=eVe2=e 

< yc,f,t+ Y yc,e,t+ Y yc,e,t+ Y yc,ei,e2,t + 

eeV+(f) eeVc'if) (ei,e2)eP<r-,(/) 

Y yc,eue2,t+ Y yc,eue2,t for a G A / £ (12) 

(ei,e2)GP+2(/) (ei-e2)eP-2(/) 

• Ordering implications 

Y x1 < |y(A)| - 1 for ah cycles A G G^'^^'^ (13) 
aey(A) 

Constraints (8), (9), and (10) represent the flow constraints for the G2SC network. Con- 
straints (11) and (12) link the action effects and prevail conditions with the corresponding 
flows, and constraint 13 ensures that the actions can be linearized into some feasible order- 
ing. 

3.5 State Change Path (PathSC) Formulation 

There are several ways to generalize the network representation of the CISC formulation 
and loosen the interaction between the networks. The G/cSC formulation presented one 
generalization that allows up to k transitions in each state variable per plan period. Since 
it uses exponentially many variables another way to generalize the network representation 
of the CISC formulation is by requiring that each value can be true at most once per plan 
period. To illustrate this idea we consider our logistics example again, but we now use 
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Load at loci 
Drive locl^locl 




Figure 10: Logistics example represented by network flow problems that allow a path of 
value transitions per plan period such that each value can be true at most once. 



a network representation that allows a path of transitions per plan period as depicted in 
Figure 10. 

Recall that the solution to the logistics example consists of three actions: first load the 
package at location 1, then drive the truck from location 1 to location 2, and last unload 
the package at location 2. Clearly, this solution would not be allowed within a single plan 
period under Graphplan-style parallelism. Moreover, it would also not be allowed within 
a single period in the GISC formulation. The reason for this is that the number of value 
changes in the package-location state variable is two. First, it changes from pack-at-locl to 
pack-in-truck, and then it changes from pack-in-truck to pack-at-loc2. As before, however, 
there does exists an ordering of the three actions that is feasible. The key idea behind this 
example is to show that we can allow multiple value changes in a single period. If we limit 
the value changes in a state variable to simple paths, that is, in one period each value is 
visited at most once, then we can still use implied precedences to determine the ordering 
restrictions. 

3.5.1 State Change Network 

In this formulation each value can be true at most once in each plan period, hence the 
number of value transitions for each plan period is limited to kc where kc = \Vc\ — 1 ioi each 
c G C. In the PathSC network, nodes appear in layers and correspond to the values of the 
state variable. However, each layer now consists of twice as many nodes. If we set up an 
IP encoding with a maximum number of plan periods T then there will be T layers. Arcs 
within a layer correspond to transitions or to value persistence, and arcs between layers 
ensure that all plan periods are connected to each other. 

Figure 11 displays a network corresponding to the state variable c with domain Vc = 
{f,g,h} that allows multiple transitions per plan period. The arcs pointing rightwards 
correspond to the persistence of a value, while the arcs pointing leftwards correspond to the 
value changes. If more than one plan period is needed the curved arcs pointing rightwards 
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link the layers between two consecutive plan periods. Note that with unit capacity on the 
arcs, any path in the network can visit each node at most once. 

PathSC network 




Period t 

Figure 11: Path state change (PathSC) network. 



3.5.2 Variables 

We now have action execution variables and arc flow variables (as deflned in the previous 
formulations), and linking variables that connect the networks between two consecutive 
time periods. These variables are defined as follows: 

• Zcj^t G {0, 1}, for c G C, / G Fc, < t < T; Zcj,t is equal to 1 if the value / of state 
variable c is the end value at period t, and otherwise. 

3.5.3 Constraints 

As in the previous formulations, we have state change flow constraints, action implication 
constraints, and ordering implication constraints. The main difference is the underlying 
network. The PathSC integer programming formulation is given as follows: 

• State change flows for all c G C, / G 14 

z = [ ^ f ^ -^oW (14) 
'^'^'^ \ otherwise. ^ ' 

yc,e,t + Zc,f,t-1 = ycj,t (15) 

yc,f,t = Yl y^'^'^ + ^'^'f'i for 1 < i < r - 1 (16) 

z,j,T = 1 iifes^c] (17) 

• Action implications for allcGC, 1 < t <T 

Y = yc,e,t for eeEc (18) 

x'i < %j,t ioTfeV,- (19) 
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• Ordering implications 



xl < |y(A)| -1 for all cycles A eG^'^^^' 



(20) 



aeV{A) 



Constraints (14)-(17) are the network flow constraints. For each node, except for the initial 
and goal state nodes, they ensure a balance of flow (i.e. flow-in must equal flow-out). The 
initial state node has a supply of one unit of flow and the goal state node has a demand of 
one unit of flow, which are given by constraints (14) and (17) respectively. The interactions 
that actions impose upon different state variables are represented by the action implication 
constraints (18) and (19), which have been discussed earlier. 

The implied precedence graph for this formulation is given by G^^^'^ = {yp^'^'^ ^^prec y 
It has an extra set of arcs to incorporate the implied precedences that are introduced when 
two actions imply a state change in the same class c € C The nodes a € V^^'^'^ again 
correspond to actions, and there is an arc (a, 6) G E^'^^^ if action a has to be executed 
before action h in the same time period, or if h has to be executed after a. More specifically, 
we have 



As before, the ordering implication constraints (20) ensure that the actions in the solu- 
tion plan can be linearized into a feasible ordering. 

4. Branch-and-Cut Algorithm 

IP problems are usually solved with an LP-based branch-and-bound algorithm. The basic 
structure of this technique involves a binary enumeration tree in which branches are pruned 

according to bounds provided by the LP relaxation. The root node in the enumeration tree 
represents the LP relaxation of the original IP problem and each other node represents a 
subproblem that has the same objective function and constraints as the root node except 
for some additional bound constraints. Most IP solvers use an LP-based branch-and-bound 
algorithm in combination with various preprocessing and probing techniques. In the last 
few years there has been significant improvement in the performance of these solvers (Bixby, 



In an LP-based branch-and-bound algorithm, the LP relaxation of the original IP prob- 
lem (the solution to the root node) will rarely be integer. When some integer variable x 
has a fractional solution v we branch to create two new subproblems, such that the bound 
constraint x < [v\ is added to the left-child node, and x > \v~\ is added to the right-child 
node. This branching process is carried out recursively to expand those subproblems whose 
solution remains fractional. Eventually, after enough bounds are placed on the variables, an 
integer solution is found. The value of the best integer solution found so far, Z*, is referred 
to as the incumbent and is used for pruning. 






{a,b)eAxA,ceC,feVc,eeES,e'eE^: 

eev+{f)/\e'eVc-{f) 



2002). 
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In a minimization problem, branches emanating from nodes whose solution value Zlp is 
greater than the current incumbent, Z*, can never give rise to a better integer solution as 
each child node has a smaller feasible region than its parent. Hence, we can safely eliminate 
such nodes from further consideration and prune them. Nodes whose feasible region have 
been reduced to the empty set, because too many bounds are placed on the variables, can 
be pruned as well. 

When solving an IP problem with an LP-based branch-and-bound algorithm we must 
consider the following two decisions. If several integer variables have a fractional solution, 
which variable should we branch on next, and if the branch we are currently working on 
is pruned, which subproblem should we solve next? Basic rules include use the "most 
fractional variable" rule for branching variable selection and the "best objective value" rule 
for node selection. 

For our formulations a standard LP-based branch-and-bound algorithm approach is very 
ineffective due to the large number (potentially exponentially many) ordering implication 
constraints in the GISC, G2SC, and PathSC formulations. While it is possible to reduce 
the number of constraints by introducing additional variables (Martin, 1991), the resulting 
formulations would still be intractable for all but the smallest problem instances. Therefore, 
we solve the IP formulations with a so-called branch- and- cut algorithm, which considers the 
ordering implication constraints implicitly. A branch-and-cut algorithm is a branch-and- 
bound algorithm in which certain constraints are generated dynamically throughout the 
branch-and-bound tree. A flowchart of our branch-and-cut algorithm is given in Figure 12. 

If, after solving the LP relaxation, we are unable to prune the node on the basis of the 
LP solution, the branch-and-cut algorithm tries to find a violated cut, that is, a constraint 
that is valid but not satisfied by the current solution. This is also known as the separation 
problem. If one or more violated cuts are found, the constraints are added to the formulation 
and the LP is solved again. If none are found, the algorithm creates a branch in the 
enumeration tree (if the solution to the current subproblem is fractional) or generates a 
feasible solution (if the solution to the current subproblem is integral). 

The basic idea of branch-and-cut is to leave out constraints from the LP relaxation 
of which there are too many to handle efficiently, and add them to the formulation only 
when they become binding at the solution to the current LP. Branch-and-cut algorithms 
have successfully been applied in solving hard large-scale optimization problems in a wide 
variety of applications including scheduling, routing, graph partitioning, network design, 
and facility location problems (Caprara &; Fischetti, 1997). 

In our branch-and-cut algorithm we can stop as soon as we find the first feasible solution, 
or we can implicitly enumerate all nodes (through pruning) and find the optimal solution 
for a given objective function. Note that our formulations can only be used to find bounded 
length optimal plans. That is, find the optimal plan given a plan period (i.e. a bounded 
length). In our experimental results, however, we focus on finding feasible solutions. 

4.1 Constraint Generation 

At any point during runtime that the cut generator is called we have a solution to the 
current LP problem, which consists of the LP relaxation of the original IP problem plus 
any added bound constraints and added cuts. In our implementation of the branch-and-cut 
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Cut generator 




Branching 



Figure 12: Flowchart of our branch-and-cut algorithm. For finding any feasible solution (i.e. 

optimize = no) the algorithm stops as soon as the first feasible integer solution 
is found. When searching for the optimal solution (i.e. optimize = yes) for the 
given formulation we continue until no open nodes are left. 
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algorithm, we start with an LP relaxation in which the ordering implication constraints are 
omitted. So given a solution to the current LP relaxation, which could be fractional, the 
separation problem is to determine whether the solution violates one of the omitted ordering 
implication constraints. If so, we identify the violated ordering implication constraints, add 
them to the formulation, and resolve the new problem. 

4.1.1 Cycle Identification 

In the GISC, G2SC, and PathSC formulations an ordering implication constraint is vio- 
lated if there is a cycle in the implied precedence graph. Separation problems involving 
cycles occur in numerous applications. Probably the best known of its kind is the travel- 
ing salesman problem in which subtours (i.e. cycles) are identified and subtour elimination 
constraints are added to the current LP. Our algorithm for separating cycles is based on 
the one described by Padberg and Rinaldi (1991). We are interested in finding the shortest 
cycle in the implied precedence graph, as the shortest cycle cuts off more fractional extreme 
points. The general idea behind this approach is as follows: 

1. Given a solution to the LP relaxation, determine the subgraph Gt for plan period t 
consisting of all the nodes a for which xf > 0. 

2. For all the arcs (a, b) G Gt, define the weights Wa,b '■= + x\ — 1. 

3. Determine the shortest path distance (ia,6 for all pairs ((a, 6) G Gt) based on arc 
weights b := 1 — (for example, using the Floyd- Warshall all-pairs shortest path 
algorithm) . 

4. If da, 6 — W},,a < for some arc (a, 6) € Gt, there exists a violated cycle constraint. 

While the general principles behind branch-and-cut algorithms are rather straightfor- 
ward, there are a number of algorithmic and implementation issues that may have a signif- 
icant impact on overall performance. At the heart of these issues is the trade-off between 
computation time spent at each node in the enumeration tree and the number of nodes 
that are explored. One issue, for example, is to decide when to generate violated cuts. 
Another issue is which of the generated cuts (if any) should be added to the LP relaxation, 
and whether and when to delete constraints that were added to the LP before. In our 
implementation, we have only addressed these issues in a straightforward manner: cuts are 
generated at every node in the enumeration tree, the first cut found by the algorithm is 
added, and constraints are never deleted from the LP relaxation. However, given the po- 
tential of more advanced strategies that has been observed in other applications, we believe 
there still may be considerable room for improvement. 

4.1.2 Example 

In this section we will show the workings of our branch-and-cut algorithm on the GISC 
formulation using a small hypothetical example involving two state variables ci and C2, five 
actions A\, A2, A3, A4, and A5, and one plan period. In particular we will show how the 
cycle detection procedure works and how an ordering implication constraint is generated. 
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Figure 13 depicts a solution to the current LP of the planning problem. For state variable 
ci we have that actions Al and A2 have a prevail condition on g, AA has a prevail condition 
on h, and action A2i has an effect that changes g into h. Likewise, for state variable C2 we 
have that action ^4 has an effect that changes g into /, action Ah changes g into h, and 
action Al has a prevail condition on /. Note that the given solution is fractional. Therefore 



some of the action variables have fractional values. In particular, we have x 



Al 



„A4 



0. 



x-^^ = 0.2, and x"^^ = x^^ = 1. In other words, actions A2 and A3 are fully executed while 
actions Al, A4 and A5 are only fractionally executed. Clearly, in automated planning the 
fractional execution of an action has no meaning whatsoever, but it is very common that 
the LP relaxation of an IP formulation gives a fractional solution. We simply try to show 
that we can find a violated cut even when we have a fractional solution. Also, note that the 
actions A4 and A5 have interfering effects in C2. While this would generally be infeasible, 
the actions are executed only fractionally, so this is actually a feasible solution to the LP 
relaxation of the IP formulation. 

State variable 2 




Al =A4 = 0.S,A2=A3 = l 




® ^ -<E) 

Al =^4 = 0.8,^5 = 0.2 



Figure 13: Solution to a small hypothetical planning example. The solution to the current 
LP has flows over the indicated paths and executes actions Al, A2, A3, A4, and 
A5. 



In order to determine whether the actions can be linearized into a feasible order- 
ing we first create the implied precedence graph G'^'^^'^ = (^yp'^''^ ^P'ec^_ where we have 
yprec ^ {Al, A2, Ai, A4, Ab) and EP''^'' = {{Al,A3),iA2,A3),{A3,AA),{AA,Al)}. The or- 
dering {A1,A3), for example, is established by the effects of these actions in state variable 
ci. Al has a prevail condition g in ci while ^3 changes g to h in ci, which implies that 
Al must be executed before ^3. The other orderings are established in a similar way. The 
complete implied precedence graph for this example is given in Figure 14. 

The cycle detection algorithm gets the implied precedence graph and the solution to the 
current LP as input. Weights for each arc (a, b) G ^P''ec g^j.g determined by the values of the 
action variables in the current solution. We have the LP solution that is given in Figure 
13, so in this example we have wai,A3 = 'WA3,A4 = 0.8, wa2,A3 = 1, and wa4,ai = 0.6. The 
length of the shortest path from Al to AA using weights Wa,b is equal to 0.4 (0.2-1-0.2). Hence, 
we have d^i^Ai = 0.4 and t«A4,Ai = 0.6. Since d^i^Ai ^ wa4_,ai < 0, wc have a violated cycle 
(i.e. violated ordering implication) that includes all actions that are on the shortest path 
from Al to AA (i.e. Al, A3, and ^14, which can be retrieved by the shortest path algorithm). 
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This generates the following ordering implication constraint Xi +XY +xf < 2, which will 
be added to the current LP. Note that this ordering constraint is violated by the current 
LP solution, as + + Xi^ = 0.8 + 1 + 0.8 = 2.6. Once the constraint is added to the 
LP, the next solution will select a set of actions that does not violate the newly added cut. 
This procedure continues until no cuts are violated and the solution is integer. 



Implied precedence graph 




Figure 14: Implied precedence graph for this example, where the labels show (1^0,61 ^a.b)- 



5. Experimental Results 

The described formulations are based on two key ideas. The first idea is to decompose the 
planning problem into several loosely coupled components and represent these components 
by an appropriately defined network. The second idea is to reduce the number of plan 
periods by adopting diff'crcnt notions of parallelism and use a branch-and-cut algorithm 
to dynamically add constraints to the formulation in order to deal with the exponentially 
many action ordering constraints in an efficient manner. 

To evaluate the tradeoffs of allowing more fiexible network representations we compare 
the performance of the one state change (ISC) formulation, the generalized one state change 
formulation (GISC), the generalized two state change (G2SC) formulation, and the state 
change path (PathSC) formulation. For easy reference, an overview of these formulations 
is given in Figure 15. 

In our experiments we focus on finding feasible solutions. Note, however, that our 
formulations can be used to do bounded length optimal planning. That is, given a plan 
period (i.e. a bounded length), find the optimal solution. 

5.1 Experimental Setup 

To compare and analyze our formulations we use the STRIPS domains from the second 
and third international planning competitions (IPC2 and IPC3 respectively). That is, 
Blocksworld, Logistics, Miconic, Freecell from IPC2 and Depots, Drivcrlog, Zenotravel, 
Rovers, Satellite, and Freecell from IPCS. We do not compare our formulations on the 
STRIPS domains from IPC4 and IPC5 mainly because of a peripheral limitation of the 
current implementation of the G2SC and PathSC formulations. In particular, the G2SC 
formulation cannot handle operators that change a state variable from an undefined value 
to a defined value, and the PathSC formulation cannot handle such operators if the domain 
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ISC 

Each state variable can change or 
prevail a value at most once per 
plan period. 

Q Q 
® ® 


CISC 

Each state variable can change (and 
prevail a value before and after each 
change) at most once per plan period. 

® ® 


G2SC 

Each state variable can change (and 
prevail a value before and after each 
change) at most twice per plan 
period. Cyclic changes {f,g,f) are 
allowed only if / is not the 
prevail condition of some action 

(h) ® ^—KA) 


PathSC 

The state variable can change any 
number of times, but each value 
can be true at most once per plan 
period. 



Figure 15: Overview of the ISC, CISC, G2SC, and PathSC formulations. 



size of the state variable is larger than two. Because of these limitations we could not test 
the G2SC formulation on the Miconic, Satellite and Rovers domains, and we could not test 
the PathSC formulation on the Satellite domain. 

In order to setup our formulations wc translate a STRIPS planning problem into a 
multi- valued state description using the translator of the Fast Downward planner (Helmert, 
2006). Each formulation uses its own network representation and starts by setting the 
number of plan periods T equal to one. We try to solve this initial formulation and if no 
plan is found, T is increased by one, and then try to solve this new formulation. Hence, 
the IP formulation is solved repeatedly until the first feasible plan is found or a 30 minute 
time limit (the same time limit that is used in the international planning competitions) is 
reached. We use CPLEX 10.0 (ILOG Inc., 2002), a commercial LP/IP solver, for solving 
the IP formulations on a 2.67GHz Linux machine with 1GB of memory. 

We set up our experiments as follows. First, in Section 5.2 wc provide a brief overview of 
our main results by looking at aggregated results from IPC2 and IPCS. Second, in Section 
5.3, we give a more detailed analysis on our loosely coupled encodings for planning and 
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focus on the tradeoffs of reducing the number of plan periods to solve a planning problem 
versus the increased difficulty in merging the solutions to the different components. Third, 
in Section 5.4 we briefly touch upon how different state variable representations of the same 
planning problem can influence performance. 

5.2 Results Overview 

In this general overview we compare our formulations to the following planning systems: 
Optiplan (van den Briel & Kambhampati, 2005), SATPLAN04 (Kautz, 2004), SATPLAN06 
(Kautz & Selman, 2006), and Satplanner (Rintanen et al., 2006)^. 

Optiplan is an integer programming based planner that participated in the optimal track 
of the fourth international planning competition^. Like our formulations, Optiplan models 
state transitions but it does not use a factored representation of the planning domain. 
In particular, Optiplan represents state transitions in the atoms of the planning domain, 
whereas our formulations use multi-valued state variables. Apart from this, Optiplan is 
very similar to the ISC formulation as they both adopt the Graphplan-style parallelism. 

SATPLAN04, SATPLAN06, and Satplanner are satisfiability based planners. SAT- 
PLAN04 and SATPLAN06 are versions of the well known system SATPLAN (Kautz & 
Selman, 1992), which has a long track record in the international planning competitions. 
Satplanner has not received that much attention, but is among the state-of-the-art in plan- 
ning as satisfiability. Like our formulations Satplanner generalizes the Graphplan-style 
parallelism to improve planning efficiency. 

The main results are summarized by Figure 16. It displays aggregate results from IPC2 
and IPCS, where the number of instances solved (y-axis) is drawn as a function of log time 
(x-axis). We must note that the graph with the IPC2 results favors the PathSC formulation 
over all other planners. However, as we will see in Section 5.3, this is mainly a reflection 
of its exceptional performance in the Miconic domain rather than its overall performance 
in IPC2. Morever, the graph with the IPCS results does not include the Satellite domain. 
We decided to remove this domain, because wc could not run it on the public versions of 
SATPLAN04 and SATPLAN06 nor the G2SC and PathSC formulations. While the results 
in Figure 16 provide a rather coarse overview, they sum up the following main findings. 

• Factored planning using loosely coupled formulations helps improve performance. Note 
that all integer programming formulations that use factored representations, that is 
ISC, CISC, G2SC, and PathSC (except the G2SC formulation which could not be 

3. We note that that SATPLAN04, SATPLAN06, Optiplan, and the ISC formulation are "step-optimal" 
while the GISC, G2SC, and PathSC formulations are not. There is, however, considerable controversy 

in the planning community as to whether the stcp-optimality guaranteed by Graphplan-style planners 
has any connection to plan quality metrics that users would be interested in. We refer the reader to 
Kambhampati (2006) for a longer discussion of this issue both by us and several prominent researchers 
in the planning community. Given this background, we believe it is quite reasonable to compare our for- 
mulations to step-optimal approaches, especially since our main aim here is to show that IP formulations 
have come a long way and that they can be made competitive with respect to SAT-based encodings. 
This in turn makes it worthwhile to consider exploiting other features of IP formulations, such as their 
amenability to a variety of optimization objectives as we have done in our recent work (van den Briel 
et al., 2007). 

4. A list of participating planners and their results is available at http://ipc04.icaps-conference.org/ 
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tested on all domains), are able to solve more problem instances in a given amount 
of time than Optiplan, which does not use a factored representation. Especially, the 
difference between ISC and Optiplan is remarkable as they both adopt the Graphplan- 
stylc parallelism. In Section 5.3, however, we will see that Optiplan does perform well 
in domains that are either serial by design or have a significant serial component. 

• Decreasing the encoding size by relaxing the Graphplan- style parallelism helps improve 
performance. This is not too surprising, Dimopoulos et al. (1997) already note that a 
reduction in the number of plan periods helps improve planning performance. How- 
ever, this does not always hold because of the tradeoff between reducing the number 
plan periods versus the increased difficulty in merging the solutions to the different 
components. In Section 5.3 we will see that different relaxations of Graphplan-style 
parallelism lead to different results. For example, the PathSC formulation shows supe- 
rior performance in Miconic and Driverlog, but does poorly in Blocksworld, Freecell, 
and Zenotravel. Likewise, the G2SC formulation does well in Freecell, but it does not 
seem to excel in any other domain. 

• Planning as integer programming shows new potential. The conventional wisdom in 
the planning community has been that planning as integer programming cannot com- 
pete with planning as satisfiability or constraint satisfaction. In Figure 16, however, 
we see that the ISC, CISC and PathSC formulation can compete quite well with 
SATPLAN04. While SATPLAN04 is not state-of-the-art in planning as satisfiabil- 
ity anymore, it does show that planning as integer programming has come a long 
way. The fact that IP is competitive allows us to exploit its other virtues such as 
optimization (Do et al., 2007; Benton et al., 2007; van den Briel et al, 2007). 




Figure 16: Aggregate results of the second and third international planning competitions. 



5.3 Comparing Loosely Coupled Formulations for Planning 

In this section we compare our IP formulations and try to evaluate the benefits of allowing 
more flexible network representations. Specifically, we are interested in the effects of reduc- 
ing the number of plan periods required to solve the planning problem versus dealing with 
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merging solutions to the different components. Reducing the number of plan periods can 
lead to smaller encodings, which can lead to improved performance. However, it also makes 
the merging of the loosely coupled components harder, which could worsen performance. 

In order to compare our formulations we will analyze the following two things. First, we 
examine the performance of our formulations by comparing their solution times on problem 
instances from IPC2 and IPCS. In this comparison we will include results from Optiplan 
as it gives us an idea of the differences between a formulation based on Graphplan and 
formulations based on loosely coupled components. Moreover, it will also show us the 
improvements in IP based approaches for planning. Second, we examine the number of 
plan periods that each formulation needs to solve each problem instance. Also, we will look 
at the tradeoffs between reducing the number of plan periods and the increased difficulty 
in merging the solutions of the loosely coupled components. In this comparison we will 
include results from Satplanner because, just like our formulations, it adopts a generalized 
notion of the Graphplan-style parallelism. 

We use the following figures and table. Figure 17 shows the total solution time (y-axis) 
needed to solve the problem instances (x-axis). Figure 18 shows the number of plan periods 
(y-axis) to solve the problem instances (x-axis), and Table 1 shows the number of ordering 
constraints that were added during the solution process, which can be seen as an indicator 
of the merging effort. The selected problem instances in Table 1 represent the five largest 
instances that could be solved by all of our formulations (in some domains, however, not 
all formulations could solve at least five problem instances). 

The label GPsteps in Figure 18 represents the number of plan steps that SATPLAN06, 
a state-of-the-art Graphplan-based planner, would use. In the Satellite domain, however, 
we use the results from the ISC formulation as we were unable to run the public version 
of SATPLAN06 in this domain. We like to point out that Figure 18 is not intended to 
favor one formulation over the other, it simply shows that it is possible to generate encod- 
ings for automated planning that use drastically fewer plan periods than Graphplan-based 
encodings. 

5.3.1 Results: Planning Performance 

Blocksworld is the only domain in which Optiplan solves more problems than our formu- 
lations. In Zenotravel and Satellite, Optiplan is generally outperformed with respect to 
solution time, and in Rovers and Freecell, Optiplan is generally outperformed with respect 
to the number of problems solved. As for the other IP formulations, the GISC provides 
the overall best performance and the performance of the PathSC formulation is somewhat 
irregular. For example, in Miconic, Driver log and Rovers the PathSC formulation does very 
well, but in Depots and Freecell it does rather poorly. 

In the Logistics domain all formulations that generalize the Graphplan-style parallelism 
(i.e. GISC, G2SC, and PathSC) scale better than the ISC formulation and Optiplan, which 
adopt the Graphplan-style parallelism. Among CISC, G2SC, and PathSC formulations 
there is no clear best performer, but in the larger Logistics problems the GISC formulation 
seems to do slightly better. The Logistics domain provides a great example of the tradeoff 
between flexibility and merging. By allowing more actions to be executed in each plan 
period, generally shorter plans (in terms of number of plan periods) are needed to solve 
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the planning problem (see Figure 18), but at the same time merging the solutions to the 
individual components will be harder as one has to respect more ordering constraints (see 
Table 1). 

Optiplan versus ISC. If we compare the ISC formulation with Optiplan, we note 
that Optiplan fares well in domains that are either serial by design (Blocksworld) or in 
domains that have a significant serial aspect (Depots) . We think that Optiplan's advantage 
over the ISC formulation in these domains is due to the following two possibilities. First, 
our intuition is that in serial domains the reachability and relevance analysis in Graphplan is 
stronger in detecting infeasible action choices (due to mutex propagation) than the network 
flow restrictions in the ISC formulation. Second, it appears that the state variables in these 
domains are more tightly coupled (i.e. the actions have more effects, thus transitions in one 
state variable are coupled with several transitions in other state variables) than in most 
other domains, which may negatively affect the performance of the ISC formulation. 

ISC versus GISC. When comparing the ISC formulation with the CISC formulation 
we can see that in all domains, except in Blocksworld and Miconic, the CISC formulation 

solves at least as many problems as the ISC formulation. The results in Blocksworld are 
not too surprising and can be attributed to semantics of this domain. Each operator in 
Blocksworld requires one state change in the state variable of the arm (stack and putdown 
change the status of the arm to arm — empty, and unstack and pickup change the status 
of the arm to holding — x where x is the block being lifted). Since, the ISC and the 
CISC formulations both allow at most one state change in each state variable, there is no 
possibility for the GISC formulation to allow more than one action to be executed in the 
same plan period. Given this, one may think that the ISC and CISC formulations should 
solve at least the same number of problems, but in this case the prevail constraints (5) of 
the ISC formulation are stronger than the prevail constraints (6) of the GISC formulation. 
That is, the right-hand side of (6) subsumes (i.e. allows for a larger feasible region in the LP 
relaxation) than the right-hand side of (5). In Figure 17 we can see this slight advantage 
of ISC over GISC in the Blocksworld domain. 

The results in the Miconic domain are, on the other hand, not very intuitive. We 
would have expected the GISC formulation to solve at least as many problems as the ISC 
formulation, but this did not turn out to be the case. One thing we noticed is that in this 
domain the GISC formulation required a lot more time to determine that there is no plan 
for a given number of plan periods. 

GISC versus G2SC and PathSC. Table 1 only shows the five largest problems in 
each domain that were solved by the formulations, yet it is representative for the whole set 
of problems. The table indicates that when Graphplan-style parallelism is generalized, more 
ordering constraints are needed to ensure a feasible plan. On average, the G2SC formulation 
includes more ordering constraints than the CISC formulation, and the PathSC formulation 
in its turn includes more ordering constraints than the G2SC formulation. The performance 
of these formulations as shown by Figure 17 varies per planning domain. The PathSC 
formulation does well in Miconic and Driverlog, the G2SC formulation does well in Freecell, 
and the CISC does well in Zenotravel. Because of these performance differences, we believe 
that the ideal amount of flexibility in the generalization of Graphplan-style parallelism is 
different for each planning domain. 
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Figure 17: Solution times in the planning domains of the second and third international 
planning competition. 
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that the formulation could not solve the problem instance within 30 minutes. 
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5.3.2 Results: Number of Plan Periods 

In Figure 18, we see that in all domains the flexible network representation of the GISC 
formulation is slightly more general than the 1-linearization semantics that is used by Sat- 
planner. That is, the number of plan periods required by the GISC formulation is always 
less than or equal to the number of plan periods used by Satplanner. Moreover, the flexible 
network representation of the G2SC and PathSC formulations are both more general than 
the one used by the GISC formulation. One may think that the network representation 
of the PathSC formulation should provide the most general interpretation of action paral- 
lelism, but since the G2SC network representations allows some values to change back to 
their original value in the same plan period this is not always the case. 

In the domains of Logistics, Freecell, Miconic, and Driverlog, the PathSC never required 
more than two plan periods to solve the problem instances. For the Miconic domain this 
is very easy to understand. In Miconic there is an elevator that needs to bring travelers 
from one floor to another. The state variables representation of this domain has one state 
variable for the elevator and two for each traveler (one to represent whether the traveler has 
boarded the elevator and one to represent whether the traveler has been serviced). Clearly, 
one can devise a plan such that each value of the state variable is visited at most twice. The 
elevator simply could visit all floors and pickup all the travelers, and then visit all floors 
again to bring them to their destination floor. 

5.4 Comparing DifTerent State Variable Representations 

An interesting question is to find out whether different state variable representations lead to 
different performance benefits. In our loosely coupled formulations we have components that 
represent multi-valued state variables. However, the idea of modeling value transitions as 
flows in an appropriately defined network can be applied to any binary or multi- valued state 
variable representation. In this section we concentrate on the efficiency tradeoffs between 
binary and multi-valued state descriptions. As there are generally fewer multi-valued state 
variables than binary atoms needed to describe a planning problem, we can expect our 
formulations to be more compact when they use a multi-valued state description. For this 
comparison we only concentrate on the CISC formulation as it showed the overall best 
performance among our formulations. In our recent work (van den Briel, Kambhampati, & 
Vossen, 2007) we analyze different state variable representations in more detail. 

Table 2 compares the encoding size for the GISC formulation on a set of problems 
using either a binary or multi-valued state description. The table clearly shows that the 
encoding size becomes significantly smaller (both before and after CPLEX presolve) when 
a multi-valued state description is used. The encoding size before presolve gives an idea of 
the impact of using a more compact multi-valued state description, whereas the encoding 
size after presolve shows how much preprocessing can be done by removing redundancies 
and substituting out variables. 

Figure 19 shows the total solution time (y-axis) needed to solve the problem instances 
(x-axis). Since we did not make any changes to the GISC formulation, the performance 
differences are the result of using different state descriptions. In several domains the multi- 
valued state description shows a clear advantage over the binary state description when 
using the CISC formulation, but there are also domains in which the multi-valued state 



250 



Loosely Coupled Formulations for Automated Planning 



description does not provide too much of an advantage. In general, however, the GISC 
formulation using a multi-valued state description leads to the same or better performance 
than using a binary state description. In all our tests, we encountered only one problem 
instance (Rovers pfilclO) in which the binary state description notably outperformed the 
multi-valued state description. 

6. Related Work 

There are only few integer programming-based planning systems. Bylander (1997) considers 
an IP formulation based on converting the propositional representation given by Satplan 
(Kautz & Selman, 1992) to an IP formulation with variables that take the value 1 if a 
certain proposition is true, and otherwise. The LP relaxation of this formulation is 
used as a heuristic for partial order planning, but tends to be rather time-consuming. A 
different IP formulation is given by Vossen et al. (1999). They consider an IP formulation 
in which the original propositional variables are replaced by state change variables. State 
change variables take the value 1 if a certain proposition is added, deleted, or persisted, 
and otherwise. Vossen et al. show that the formulation based on state change variables 
outperforms a formulation based on converting the propositional representation. Van den 
Briel and Kambhampati (2005) extend the work by Vossen et al. by incorporating some of 
the improvements described by Dimopoulos (2001). Other integer programming approaches 
for planning rely on domain-specific knowledge (Bockmayr & Dimopoulos, 1998, 1999) or 
explore non-classical planning problems (Dimopoulos & Gerevini, 2002; Kautz k, Walser, 
1999). 

In our formulations we model the transitions of each state variable as a separate flow 
problem, with the individual problems being connected through action constraints. The 
Graphplan planner introduced the idea of viewing planning as a network flow problem, but 
it did not decompose the domain into several loosely coupled components. The encodings 
that we described are related to the loosely-coupled modular planning architecture by Sri- 
vastava, Kambhampati, and Do (2001), as well as factored planning approaches by Amir 
and Engelhardt (2003), and Brafman and Domshlak (2006). The work by Brafman and 
Domshlak, for example, proposes setting up a separate CSP problem for handling each fac- 
tor. These individual factor plans are then combined through a global CSP. In this way, it 
has some similarities to our work (with our individual state variable flows corresponding to 
encodings for the factor plans). Although Brafman and Domshlak do not provide empirical 
evaluation of their factored planning framework, they do provide some analysis on when 
factored planning is expected to do best. It would be interesting to adapt their minimal 
tree-width based analysis to our scenario. 

The branch-and-cut concept was introduced by Grotschel, Riinger, and Reinelt (1984) 
and Padberg and Rinaldi (1991), and has been successfully applied in the solution of 
many hard large-scale optimization problems (Caprara & Fischetti, 1997). In planning, 
approaches that use dynamic constraint generation during search include RealPlan (Srivas- 
tava et al., 2001) and LPSAT (Wolfman & Weld, 1999). 

Relaxed definitions for Graphplan-stylc parallelism have been investigated by several 
other researchers. Dimopoulos et al. (1997) were the first to point out that it is not necessary 
to require two actions to be independent in order to execute them in the same plan period. 
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Table 2: Formulation size for binary and multi-valued state description of problem instances 
from the IPC2 and IPCS in number of variables (^va), number of constraints 
(#co), and number of ordering constraints, or cuts, (#cu) that were added dy- 
namically through the solution process. 
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Figure 19: Comparing binary state descriptions with multi-valued state descriptions using 
the GISC formulation. 
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In their work they introduce the property of post-serializabiUty of a set of actions. A set 
of actions A' Q A is said to be post-scriahzable if (1) the union of their preconditions 
is consistent, (2) the union of their effects is consistent, and (3) the preconditions-effects 
graph is acychc. Where the preconditions-effects graph is a directed graph that contains 
a node for each action in A', and an arc (a, b) for a,b e A' if the preconditions of a are 
inconsistent with the effects of a. For certain planning problems Dimopoulos et al. (1997) 
show that their modifications reduce the number of plan periods and improve performance. 
Rintanen (1998) provides a constraint-based implementation of their idea and shows that 
the improvements hold over a broad range of planning domains. 

Cayrol et al. (2001) introduce the notion of authorized linearizations, which implies an 
order for the execution of two actions. In particular, an action a € A authorizes an action 
b e A implies that if a is executed before b, then the preconditions of b will not be deleted 
by a and the effects of a will not be deleted by b. The notion of authorized linearizations 
is very similar to the property of post-serializability. If we were to adopt these ideas in 
our network-based representations it would compare to the GISC network in which the 
generalized state change arcs (see Figure 7) only allows values to prevail after, but not 
before, each of the transition arcs. 

A more recent discussion on the definitions of parallel plans is given by Rintanen, Hel- 
janko and Niemela (2006). Their work introduces the idea of 3-step semantics, which says 
that it is not necessary that all parallel actions are non-interfering as long as they can be 
executed in at least one order. 3-step semantics provide a more general interpretation of 
parallel plans than the notion of authorized linearizations used by LCGP (Cayrol et al., 
2001). The current implementation of 3-step semantics in Satplanner is, however, some- 
what restricted. While the semantics allow actions to have conflicting effects, the current 
implementation of Satplanner does not. 

7. Conclusions 

This work makes two contributions: (1) it improves state of the art in modeling planning as 
integer programming, and (2) it develops novel decomposition methods for solving bounded 
length (in terms of number of plan periods) planning problems. 

We presented a series of IP formulations that represent the planning problem as a set 
of loosely coupled network flow problems, where each network flow problem corresponds 
to one of the state variables in the planning domain. We incorporated different notions 
of action parallelism in order to reduce the number of plan periods needed to find a plan 
and to improve planning efficiency. The IP formulations described in this paper have led 
to their successful use in solving partial satisfaction planning problems (Do et al., 2007). 
Moreover, they have initiated a new line of work in which integer and linear programming 
are used in heuristic state-space search for automated planning (Benton et al. , 2007; van den 
Briel et al., 2007). It would be interesting to see how our approach in the context of IP 
formulations applies to formulations based on satisfiability and constraint satisfaction. 
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